Sum rules and vertex corrections for electron-phonon interactions 



O 
O 

(N 



0. RoschQ G. Sangiovanni, and 0. Gunnarsson 
Max- Planck- Institut fur Festkorperforschung, Heisenbergstr. 1, D-70506 Stuttgart, Germany 

We derive sum rules for the phonon self-energy and the electron-phonon contribution to the 
electron self-energy of the Holstein-Hubbard model in the limit of large Coulomb interaction U. 
Their relevance for finite U is investigated using exact diagonalization and dynamical mean-field 
theory. Based on these sum rules, we study the importance of vertex corrections to the electron- 
phonon interaction in a diagrammatic approach. We show that they are crucial for a sum rule for 
the electron self-energy in the undoped system while a sum rule related to the phonon self-energy 
of doped systems is satisfied even if vertex corrections are neglected. We provide explicit results for 
the vertex function of a two-site model. 
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I. INTRODUCTION 



Recently, there has been much interest in the pos- 
sibility that electron-phonon interactions may play an 
important role for properties of cuprates, e.g., for 
superconductivity. In particular, the interest has fo- 
cused on the idea that the Coulomb interaction U might 
enhance effects of electron-phonon interactions, e.g., due 
to interactions with spin fluctuations. 4 Effects of the 
electron-phonon coupling are described by the electron- 
phonon part S op of the electron self-energy £ and the 
phonon self-energy, II. We have previously derived sum 
rules for these quantities for the t-J model, and used the 
sum rules to demonstrate that the electron-phonon in- 
teraction influences S ep and II in quite different ways for 
strongly correlated systems^i Here, we extend this work 
and derive sum rules for the related half-filled Holstein- 
Hubbard model in the limit of a large U. We obtain sum 
rules for E cp integrating either over all frequencies or 
only over frequencies in the photoemission energy range. 
The latter sum rule shows a very strong dependence on 
U , increasing by a factor of four in going from U = 
to U = co. From numerical calculations, using both ex- 
act diagonalization and dynamical mean-field theory, we 
show that the U = co result is also relevant for inter- 
mediate values of U w 3-D, where D is half the band 
width. 

In a diagrammatic many-body language, the electron- 
phonon interaction could be enhanced by U via cor- 
rections to vertex functions^ or dressing of Green's 
functions^ Huang et alA and Koch and Zeyher— stud- 
ied how U changes an effective vertex function in the 
static limit {ui = 0) and found a suppression, although it 
was concluded in Ref . |5| that the suppression is reduced 
for a large U and a small wave vector q. Often, one is 
not only interested in these special cases but in proper- 
ties that depend on integrals over oj and q containing 
vertex functions. Here, we study to what extent the sum 
rules above are fulfilled when vertex corrections are ne- 
glected. We find that the sum rule for E op , integrating 
over the photoemission energy range, is underestimated 
by a factor of four if vertex corrections are neglected. On 



the other hand, a sum rule for the phonon self-energy is 
fulfilled also without vertex corrections. This suggests 
that it can be important to include vertex corrections for 
studying properties of cuprates and other strongly corre- 
lated materials. 

The Hubbard model with electron-phonon interaction 
is introduced in Sec. [HJ In Sec. IIIII sum rules for the 
electron and phonon self-energies are derived focusing on 
the limit U — > co and in Sec. IIVI we numerically check 
their accuracy for large but finite U. These sum rules 
then form the basis for the discussion of the effects of 
vertex corrections in Sec. The results are illustrated 
in Sec.lVTJfor a two-site model. 



II. 



HUBBARD MODEL IN THE LIMIT OF 
LARGE U 



Strongly correlated electrons are often described by the 
Hubbard model 
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where Ed is the level energy, t(> 0) is the hopping integral 
between nearest-neighbor sites U is the Coulomb 

repulsion between two electrons on the same site, ct cre- 
ates an electron on site i with spin a, and — c\ a c ia . 
In addition, we introduce an electron-phonon interaction 

H ^ = 7f E ~ !)( & q + ^ q )e iq - R % (2) 



where N is the number of sites, 



i lT -I- nil > and frq 
creates a phonon with the wave vector q and energy w q 
as described by the free phonon Hamiltonian 



H ph = ^w q 6 q 6 q . 



(3) 



We assume a q-dependent on-site coupling with the 
strength (? q . The coupling to hopping integrals is ne- 
glected, which, e.g., has been found to be a good approx- 
imation for the planar oxygen (half-)breathing mode in 
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the high-T c cuprates£ The special case of a Holstein cou- 



pling is obtained by setting g^—g and - 



L> p h for all 



To describe photoemission (PES) and inverse photoe- 
mission (IPES) within the sudden approximation, we 
consider the one-electron removal (-) and addition (+) 
spectra 



A-(k,u) 



A+(k,u>) = 



E- 

mn 

x -> e 



Z 



run 



\(n\c k Jm)\ 2 S(Lo+£ n -£ m ), (4) 



\(n\cljm)\ 2 S(uj+£ m -£ n ), (5) 



where the energy w is measured relative to the chemi- 
cal potential /z, \m) and \n) are eigenstates of the grand 
canonical Hamiltonian TL = H — fi(N) with eigenenergies 
£ m /n, Z = X] m exp(-^m) is the corresponding parti- 
tion sum, and j3 — l/T. We assume that there is no 
explicit dependence on the electron spin a. From the 
total spectral density 



A(k,u) =A-{k,u)+A+(k,co), 
we obtain the one-electron Green's function 

A(k,w) 



G(k,z) = [ 



(lid 



(6) 



(7) 



which depends on the electronic wave vector k and the 
complex energy z. It is related to the electron self-energy 
E(k, z) via the Dyson equation 



G(k, z) 



1 



-£k-E(k,; 



(8) 



where e k = £<i-iI](i.j)( e ' k(Ri ~ R,) +H.c.)/JV is the bare 
electronic dispersion. We can split up G(k, z) into two 
parts, 



G(k,z) = G+(k, z)+G-(k,z), 



(9) 



where the (I) PES Green's functions G ± (k, z) are defined 
by replacing the spectral density in Eq. Q by A ± (k, uj). 
We also define corresponding self-energies by 



G ± {k,z) 



E±(k, 



(10) 



where is the integrated weight and absorbs energy- 
independent contributions to the self-energy. We con- 
sider the half-filled system and choose Ed = — U/2 to 
have explicit particle-hole symmetry Then, as detailed 
in the appendix, — * 1/2, — > ±U /2 + £k for U — ► oo. 
In this limit, one finds from inserting Eqs. © and i|l(J|) 
into Eq. JjJJ to leading order m 1/U 

Im E(k, z w ±t//2) = 2 Im S ± (k, z w ±C//2), (11) 

relating the spectral densities of the different self-energies 
for energies z pa ±{7/2. 



In the limit of large f, states with double occupancy 
can be projected out. This leads to the t-J modelS. 
as an effective low-energy model for Eq. Q if certain 
terms are assumed to be negligible*^ as will be done 
here. Since double occupancy is excluded in the t-J 
model, its electron Green's function has no contribution 
from inverse photoemission for the undoped system, i.e., 
Gt-j(k,z) = G^_j(k,z). We assume that the photoe- 
mission spectra of the Hubbard and the t-J model are 
identical (apart from a trivial energy shift ps U/2), and 
thus 



G"(k, z) = G t .j(k,z + U/2). 



(12) 



The inverse photoemission in the half-filled Hubbard 
model can be related to the photoemission in the un- 
doped t-J model because of particle-hole symmetry, 



G+(k, 



G t -j(k AF -k,U/2- z) 



(13) 



for a two-dimensional square lattice where k^F = 
(tv/cl, /a). Terms of order t <§; U are neglected in the 
energy shifts ±U/2 in Eqs. lO and JT3J. 

The electron self-energy in the t-J model is defined 
via Gt-j(k, z) = 0.5/fw — E f _j(k, z)] and it follows from 
Eas. pi ty . 1(12 )1 . and ifE ty that 



E~(k,«) = Z t .j(k,z + U/2), 
E+(k,z) = E w (k A F - k, ■ 

III. SUM RULES 



z). 



(14) 
(15) 



In this section, we discuss sum rules for the electron 
and phonon self-energies. They will allow us to ad- 
dress the importance of vertex corrections to the electron- 
phonon interaction in a diagrammatic treatment of the 
Hubbard model. 



A. Electron self-energy 

In App.IaTI we use spectral moments to derive the fol- 
lowing sum rule for the electron self-energy E(k, z) which 
gives the total weight of its spectral density integrated 
over all frequencies: 

1 f°° 

- / dw Im Efj (k, uj ^0+) = U 2 (n t _ a ) (1 - (n M )) 
- q * 

+ ^ E \9«\ 2 (K + & -q| 2 ) - 4=o(^ + (16) 
q 

where we defined p qcr = J2i ^i<Te~ jqRi /V^V and 6, = 
SqV^Vv^- CT i s t ne electron spin for which the 
self-energy is calculated but in our case the results do not 
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depend on it. The sum rule in Eq. (|16fl is valid for any U 
and interestingly, it is independent of the electronic wave 
vector k. For a Holstein coupling, it simplifies to 

1 f°° 

- / ( LImE ff (k, W - ! + ) = [/ 2 (n 1 J(l-( n! _ ff )) 

K J-aa 

+2Ug((b l +b\)n l _ a )+g 2 ((b l +bl) 2 }. (17) 

In the derivation of Eqs. I|l(jfl and l|17l) . we have assumed 
translation invariance so the expectation values on the 
right hand side of the equations do not depend on the 
site i at which they are evaluated. 

We now focus on the half-filled system where the mean 
occupation per site and spin (rii_ CT ) = 1/2 and consider 
the limit U — > oo. Because of the complete suppression of 
double occupancies and the specific form of the electron- 
phonon coupling in Eq. (J5J, there are no phonons ex- 
cited in the ground state and expectation values involving 
phonon operators in Eqs. (|16fl and (|17fl greatly simplify, 
e.g., U{(h + b\)m- a ) -» and ((&; + b\) 2 } -> 1. For the 
electron-phonon contribution S cp (k, z) to the electron 
self-energy, i.e. the difference between the self-energies 
for systems with and without electron-phonon coupling, 
one then obtains the following sum rule: 



9a\ 



1 f°° 1 
lim - / duj Im S op (k, uj - i0+) = — V 

(18) 

In the special case of a Holstein coupling, g = g. 

As described in App. IjA 2(1 . corresponding sum rules 
for the electron-phonon contributions to the (I)PES self- 
energies S ± (k, z) can be derived analogously. One finds 



lim - / duj Im S ± ' ep (k,tj - i0+) = 

(7^oo 7T J_ (x> 



(19) 



It follows from Eqs. and JT^J that Eq. (0 directly 
translates into a sum rule for the electron-phonon con- 
tribution to the electron self-energy in the undoped t-J 
model. This result was obtained already in Ref. us- 
ing spectral moments in the framework of the t-J model 
itself. 

The spectral functions of S ± (k, z) are non-zero only in 
the energy range where the (I)PES spectra Im G ± (k, w — 
iQ + )/-K are located, i.e., around uj » ±{7/2. With 
Eq. (|llf> . we can therefore derive from Eq. I|19|) also par- 
tial sum rules for the electron-phonon contribution to the 
full self-energy S(k, z) in the Hubbard model, 



lim — 



(I)PES 



duj Im E ep (k, u - i0 + ) = 2g 2 , (20) 



where we integrate over the (I) PES energy range around 
uj ps ±{7/2. An explicit choice in the limit of large U 
would be, e.g., to integrate from — oo to — [7/4 (from {7/4 
to oo) in order to fully include the energy range of the 
PES (IPES) spectrum. 
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FIG. 1: Weights obtained by integrating Im E ep (k, uj—i0 + )/n 
over the indicated frequency intervals for the half-filled Hub- 
bard and undoped t-J models. Also shown are the result for 
the Hubbard model without vertex corrections (no VC) and 
the lowest-order result for the Hubbard model with U = 
(non-int.). For the t-J model, the photoemission spectrum 
has been shifted by -U /2 and for the (7 = Hubbard model 
the photoemission and inverse photoemission spectra have 
been shifted by -{7/2 and {7/2, respectively. The results for 
{7 = refer to the k-averaged self-energy. 



Together with Eq. JTBJl, it follows from Eq. (|2U| that 
the corresponding partial sum rule for integrating over 
the central energy range u> ~ is given by 



lim - f duj Im E op (k,w - iQ+) = -3g 2 

U~> CO 7T J^q 



(21) 



The negative value can be understood as follows. For 
g = 0, Im S(k, uj — i0 + )/7r has a pole in this energy 
range with a large positive weight (w (7 2 /4). When the 
electron-phonon coupling is switched on, the strength of 
this pole is slightly reduced which shows up as a pole with 
a negative weight in the spectral function of E cp (k, z). 

The different sum rules for the electron-phonon con- 
tribution to the self-energy in the Hubbard and the t- 
J model are summarized schematically in Fig. ^ F° r 
comparison, we also show the result for non-interacting 
electrons (U = 0) to lowest order in ~g 2 . Integrating the 
k-averaged S° p n _ int over the range of the photoemission 
spectrum for a half-filled model gives the contribution 
g 2 /2, which is a factor of 4 smaller than what is ob- 
tained in the largc-{7 model (cf. Eq. (J2JJJ). Results for 
the large-{7 Hubbard model obtained by neglecting ver- 
tex corrections (no VC) are also shown in Fig. ^ and will 
be discussed in Sec. 



B. Phonon self-energy 

To lowest order in <? q , the phonon self-energy n(q, z) 
is given by \gq\ 2 x(<l, z)/N where x(q, z) is the charge- 
charge response function. Sum rules for xi% z ) therefore 
directly translate into approximate sum rules for n(q, z). 
At T = 0, 



x(q^) 



i 



Z + UJy 



(22) 



4 



where \v) is an eigenstate of H with eigenenergy oj v rel- 
ative to the ground state energy and p q = /9 qcr is the 
Fourier transform of m. For U 3> \t\ and a half-filled sys- 
tem, the ground state has exactly one electron per site to 
lowest order in t/U. Applying p q to the ground state |0), 
it then follows that p q |0) is zero to this order if we con- 
sider q ^ 0. The sum rule for |Im x(q, to + i0 + )\/ (irN) 
integrated over all frequencies is then also zero to lowest 
order in t/U: 

1 f°° 

— J dio\Im X (<i^0,LU + i0 + )\=O(t/U). (23) 

In the t-J model, an exact sum rule for the spectral 
function of %(q, z) has been derived^ 

1 f°° 

~^J2 du\hn X trj(tl,cj + iO + )\=2S(l-S), (24) 

where S is the doping. We assume as before that the 
photoemission spectra of the large-C/ Hubbard and the 
t-J model are identical. They extend over an energy 
range which is equal or smaller than the width of the 
lower Hubbard band, A = 0(t) < U. Then, Eq. JUJ 
also leads to a sum rule for x(q, z) if the integration is 
limited to \u>\ < 2A which excludes transitions from the 
lower to the upper Hubbard bands not captured by the 
t-J model: 



U/D 



1 



E 



r 2A 



x(q,uj + i0 + )\=25(l-5). 



(25) 



> / duj llm 

ttN 2 ^ J 2A 1 

As discussed in Ref. 0, these results indicate a strong 
suppression of the phonon self-energy in weakly doped 
systems with strong correlations because one would ob- 
tain unity if non-interacting electrons (U = 0) were as- 
sumed instead. 



IV. NUMERICAL RESULTS AT FINITE U 



The sum rules in Sec. IIII Al for the electron-phonon 
contribution to the electron self-energy have been derived 
in the limit U — ► oo. In practice, however, we are inter- 
ested in strongly correlated systems with large but finite 
U . In order to check the usefulness of the sum rules for 
such cases we have performed numerical calculations us- 
ing exact diagonalization (ED) and dynamical mean-field 
theory (DMFT). 

The ED calculations are done on a two-dimensional 
tilted 10-site square cluster with periodic boundary con- 
ditions. We consider a weak electron-phonon coupling 
such that it suffices to include only states with at most 
one phonon excited, thereby limiting the size of the 
phonon Hilbert space. In addition to these exact re- 
sults that are only influenced by finite size effects, we 
also consider the thermodynamic limit in the dynamical 
mean-field approximation^ which in addition allows us 




■\/{U/Df -> 

FIG. 2: Relative deviation of the total sum rule for the lo- 
cal (k-averaged) E ep from its large-!/ limit (Eq. I)18ft ') as a 
function of U using ED and DMFT with A = 0.0025. 



to consider larger electron-phonon couplings. Instead of 
studying their effects in the often assumed paramagnetic 
phased we do DMFT in the antifcrromagnetic phase 
(AF-DMFT) since antiferromagnetic correlations are be- 
lieved to be important J£ We consider a Bethe lattice 
with a semi-elliptical density of states (half bandwidth 
D). The self-consistent Anderson impurity model which 
appears in this approach is solved using exact diagonal- 
ization and a continued fraction expansion^ (employ- 
ing up to 14 discrete bath levels and allowing for up 
to 30 excited phonons). Specifically, we study the un- 
doped Hubbard-Holstein model with D = At = 1 and 
LUph = 0.025D = O.lt. In both approaches, we first calcu- 
late the (I)PES spectral densities A ± which directly lead 
to the Green's function G using Eqs. JSJ and 0. The 
inversion of the Dyson equation (Eq. JSJ) then gives the 
self-energy S whose spectral density we integrate over 
different energy ranges. Alternatively, the total sum rule 
can be obtained from calculating the ground state expec- 
tation values appearing in Eq. I|17|) . 

We first consider a weak electron-phonon coupling cor- 
responding to a dimensionless coupling constant A = 
g 2 /(uj ph D) = g 2 /(w ph 4t) = 0.0025. Using results from 
both ED and DMFT calculations, Fig.|2shows how much 
the total spectral weight of the electron-phonon contri- 
bution to the local (k-averaged) electron self-energy at fi- 
nite U deviates from the sum rule for U — > oo (Eq. ^ . 
Relative to the latter, the deviation is less than 10% for 
U as small as 3.5-D, a value often considered appropriate 
for the cuprates. This difference decreases like 1/U 2 as 
can be seen from Fig. [3 in agreement with expectations 
from a simple pcrturbational approach which we discuss 
in more detail at the end of this section. We note that 
results from ED and DMFT agree rather well. This con- 
sistency indicates that the finite size effects of the former 
and the approximations of the latter method are proba- 
bly not strongly influencing the results discussed here. 

The full U dependence of both total and partial sum 
rules is illustrated in Fig. fusing results from DMFT cal- 
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FIG. 3: U dependence of total and partial sum rules for E ep 
using DMFT with A = 0.0025. The dotted lines indicate the 
expected small-?/ behavior. 
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FIG. 4: Relative deviations of total and partial sum rules 
from their large-?/ limits as a function of A using DMFT for 
U = 5D. For the total sum rule, also results from perturba- 
tion theory (PT) for a simplified model (see text) are shown. 



deviations are comparable in size to that of the total sum 
rule (this similarity is also found when the U dependence 
is considered) and are less than 10%. In all cases, the 
deviations decrease linearly with A. Again, for the to- 
tal sum rule, the results can be quite well described by 
a result from perturbation theory for a simplified model 
which we introduce in the following. 

We consider the self-consistent Anderson impurity 
model to be solved in AF-DMFT of the Hubbard- 
Holstein model and replace the bath by a single (but 
spin-dependent) level. With the impurity level at —U/2, 
the self-consistent bath level is located at sa ±U/2 de- 
pending on the spin orientation; the hopping amplitude 
between impurity and bath is fixed to V ~ D/2. We treat 
both this hybridization and the electron-phonon interac- 
tion as perturbations of the atomic limit and find for 
large U 



2V 2 g 



U 2 (U + tu ph ) 



+ 0(U- 4 ) (26) 



and 



{(Pi 



l 



W 2 g 2 /uj ph 



U{U 



o{u- 



- w p h)(C/ + 2w ph ) 

(27) 

Using these results in Eq. (|17|l . we expect the total spec- 
tral weight of the electron-phonon contribution to the 
electron self-energy to be proportional to U~ 2 as was ob- 
served in Fig. [21 Although we have replaced the bath by 
a single level, the expressions in Eqs. (|26|l and (|27|l give 
a rather accurate description of the numerical results as 
can be seen in Fig. (J3J. 

In conclusion, we find that for typical values of U the 
relative deviations from the sum rules derived for U — > oo 
are smaller than 10%, and therefore these sum rules can 
be used semi-quantitatively also for finite U. 



culations. Already for U larger than the non- interacting 
bandwidth 2D, the sum rules clearly approach their re- 
spective large-t/ limits (Eqs . (fTHfl . (f2Tj|) . and (ED)- As 
discussed at the end of Sec. IIII Al we expect for U — ► 
that the total sum rule again approaches g 2 with the 
weight equally distributed over the PES and IPES en- 
ergy ranges and no central contribution to the spectrum 
of S ep at u> = 0. These trends are manifest in the calcu- 
lated results. 

Next, we study the dependence on the strength of the 
electron-phonon coupling for fixed U = 5D. We restrict 
ourselves to DMFT calculations where larger couplings 
are accessible because of the smaller phonon Hilbert 
space compared to ED. In Fig. 0] we plot the A depen- 
dence of the relative deviations of total and partially inte- 
grated spectral weights of the electron-phonon contribu- 
tion to the local electron self-energy from their respective 
sum rules for U —> oo (Eqs. JTHJ), (|2UJ), and For 
both the (I)PES and the central sum rules, the relative 



V. EFFECTS OF VERTEX CORRECTIONS 

We now use the sum rules introduced in Sec. 11111 to 
study the effects of vertex corrections in a diagrammatic 
calculation of the electron and phonon self-energies. We 
define the vertex function T(k,q) as the sum of all irre- 
ducible vertex diagrams connecting two electron Green's 
functions with a phonon propagator taking out one cou- 
pling constant g q explicitly (see Fig. (5J. We use the 4- 
vectors k and q as shorthand notation for the momenta 
and frequencies involved. 



A. Electron self-energy 

An important lowest-order (in |g q | 2 ) contribution to 
the electron-phonon part of the electron self-energy is 
shown in Fig. 0a,, although there are also other, more 
complicated lowest-order diagrams^ The diagram in 
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FIG. 5: Diagrammatic representation of the vertex function 
T(k, q), where k (q) stands for the incoming electron (phonon) 
momentum and frequency. The full and dashed lines represent 
electron and phonon Green's functions, respectively. 



a) 




k+q 





FIG. 6: a) Lowest-order contribution (in |<? q | 2 ) to the 
electron-phonon part of the electron self-energy E cp and b) 
the charge-charge response function \- The full and dashed 
lines represent electron and phonon Green's functions, respec- 
tively, and the circles the vertex functions (<? q )r. 



where we have also used that at half-filling, the phonon 
Green's function D{q) is not dressed in the large-C/ limit. 
Comparing this result for the diagram in Fig. with- 
out vertex corrections (no VC) with the corresponding 
exact sum rule, Eq. (|20l) . shows that this approximation 
underestimates the sum rule by a factor of four. This 
result is schematically indicated in Fig. ^ When vertex 
corrections are neglected, the diagram in Fig. has no 
contributions in the energy range u> ~ 0. 

We have elsewhere^ used the self-consistent Born ap- 
proximation to study electron-phonon interaction in the 
undoped t-J model which is closely related to the half- 
filled Hubbard model in the large-t/ limit. In this ap- 
proach, fairly good agreement with exact results is ob- 
tained although vertex corrections are neglected. Al- 
ready the lowest-order (in the electron-phonon coupling) 
diagram for the electron-phonon contribution to the elec- 
tron self-energy fulfills an exact sum rule for the total 
spectral weight. The reason for this result contrasting 
the strong violation of the sum rule in the large-?/ Hub- 
bard model when vertex corrections are neglected can be 
traced to the use of a Green's function for spinless holons 
in the self-consistent Born approximation. Its spectral 
function integrates to unity over the photoemission en- 
ergy range whereas the spin-dependent electron Green's 
function in the large-!/ Hubbard model has the spectral 
weight one half in both the photoemission and the inverse 
photoemission energy range. 



B. Phonon self-energy 

For q 7^ 0, the charge-charge response function, 
x(q, of), for the Hubbard model can be obtained from 
the diagram in Fig. |Sp: 



-*£/ 



27 



G(k + q)G(k)r{k + q,-q). (30) 



Fig. HJl is 



T, CP (k,LU 

i 

N 



£KI 



duj' 
~2T 



(28) 

G(k + q)D(q)T(k,q)T(k + q,-q), 



where q stands for a wave vector q and a frequency u>'. 
G and D are fully dressed electron and phonon Green's 
functions. 

We now neglect the vertex corrections, i.e, we put 
r = 1. After using Eq. J7J) to express the electron Green's 
function in terms of its spectral function, the u' integral 
can be performed. For the half-filled Hubbard model in 
the large-// limit, the spectral function integrated over 
the lower or the upper Hubbard band gives half an elec- 
tron per spin. As a result, we find 



lim — f did Im 

U^oo 7T 7(I)PES 



.... vc( k > w - i0+ ) = \l 



(29) 



We consider the large-?/ limit for a hole-doped Hubbard 
model. The k-averaged photoemission spectrum for a 
given spin integrates to (1 — 6)/2. As earlier, we use the 
assumption that this spectrum agrees with the spectrum 
of the t-J model and that it only extends over an energy 
range A <C U. In inverse photoemission, the probability 
of adding an electron to an unoccupied site is 5. We 
therefore assume that U is so large that the integral of 
the inverse photoemission spectrum for a given spin up to 
A is given by S. We neglect vertex corrections and replace 
T(p + q, —q) by unity. Introducing spectral functions, the 
u) integration is performed. We consider the integral of 
the spectrum of the charge-charge response function over 
\u)\ < 2 A, thereby excluding transitions between the two 
Hubbard bands. Then, 

1 f 2A 

J ^ duj l Im Xno vc(q,^ + «0 + )| 
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TABLE I: Sum rules equivalent to Eq. l-'ilt but for individual 
values of q, using ED of the t-J model on an 18-site cluster 
with two holes and J/t — 0.3. The "No VC" results were 
obtained from wp(k) and wip(k), using the second line of 
Eq. 13111 . and the exact ones result from the direct calculation 
of X- "Ratio" shows the ratio of these results. 



q/f 


(1,1) 


(2,0) 


(2,2) 


(3,1) 


(3,3) 


No VC 0.1848 


0.1927 0.2103 


0.2025 


0.2285 


Exact 


0.2100 


0.1961 


0.2191 


0.2085 


0.2212 


Ratio 


0.8804 


0.9825 


0.9597 


0.9714 


1.0330 



= t Wp ( k ) w ip ( k + q) + w p ( k + q) w ip ( k )] 

q k 

= 26(1-6), (31) 

where Wjp (k) is the integrated weight of the photoemis- 
sion spectrum for the wave vector k and wip (k) is the cor- 
responding quantity for inverse photoemission, exclud- 
ing the upper Hubbard band. For a large system, the 
q = term gives a negligible contribution in Eq. (|3T| . It 
then agrees with the sum rule in Eq. derived from 

the corresponding sum rule in the t-J model, although 
vertex corrections have been neglected. It is important, 
however, to use dressed Green's functions in calculating 
x(q,oj). Otherwise, 26 in Eq. i|31|) would have been re- 
placed by (1 +6), and there would have been a strong 
disagreement with Eq. I|25[) in the low-doping regime. 

The sum rule in Eq. i|31|) refers to an average over q. 
We next study individual values of q for a, t-J model on 
a two-dimensional tilted 18-site square cluster with two 
holes, periodic boundary conditions, and J/t — 0.3. We 
have calculated wp (k) and wip (k) , using exact diagonal- 
ization. From the second line of Eq. I)31|l . we can obtain 
sum rules for each q similar to Eq. (|31l) . The result is 
shown in the line "No VC" of Table [Q The results are 
compared with the sum rule for x*-j(q, z ) calculated in 
the t-J model ("Exact" in Table Uf) . As can be seen from 
the ratios of these results, the sum rule for individual 
values of q that can be deduced from Eq. I|31|l is also 
rather well fulfilled (typically, with a deviation of 5-10%) 
although vertex corrections are neglected. 

For the half-filled Hubbard model, there is no contri- 
bution to Im x(q, w + *0 + ) f° r M ^ U. We therefore 
focus on contributions for \lo\ ps U. If |Imx no vc(q, w + 
i0 + )|/(-7riV) is integrated over all frequencies, we obtain 
unity. The exact result in Eq. 1|23|) . however, is zero to 
lowest order in t/U. This dramatic disagreement shows 
the importance of vertex corrections in this case. 

VI. EXAMPLE: TWO-SITE MODEL 

To study the vertex corrections more explicitly, we con- 
sider a two-site Hubbard model. The electron-phonon 
interaction in Eq. @ can then be split in q = and 
q = n terms. The q = term has just the rather triv- 



ial but important effect of convoluting the spectra with 
phonon satellites, while the q = n term introduces dy- 
namics, scattering electrons between bonding and anti- 
bonding orbitals. We therefore only keep the the more 
interesting q = n term here. 

Following Huang et al.^ we can calculate the vertex 
function explicitly for the two-site model in the limit of 
U/t very large. We consider an incoming electron in the 
bonding orbital (+) with the frequency to scattered by 
the antibonding (q — 7r) phonon with the frequency lo' 
into the antibonding orbital (-) with the frequency lu + lu' 
and obtain 

r(w,+;w') = T(u + u/,—,-u/) (32) 
_ lu{lo + uj') +uj't + U 2 /4 
~ (u + t)(u + u)' -t) ' 

where various terms of higher order in t/U have been 
neglected. 

Using this result for the vertex function, we can cal- 
culate the diagram in Fig. according to Eq. i|28|l . In 
the limit of large U, we find poles with weight 2g 2 both 
atw« —U/2 and at uo « U /2. Therefore, the sum rules 
for integrating over either the PES or the IPES energy 
range, Eq. (25)) ■ are exactly fulfilled. Without vertex cor- 
rections, these sum rules are underestimated by a factor 
of four, cf. Eq. (|2"5|) . This can be understood by noting 
that for electronic energies in the range of the lower or 
upper Hubbard band, \lo\ C//2, it follows from Eq. Ij32(l 
that 

r(|u;|«^,+;a/)«2 (33) 

when the phonon frequency lo' is assumed to be small 
compared to U . Therefore, including vertex corrections 
effectively increases the weight of poles around \lu\ w U/2 
by a factor T 2 = 4. In addition, the self-energy calculated 
using vertex corrections also has poles at lu = t — w p h and 
lu = —t, the latter being a double pole. Except for a dif- 
ferent sign, to leading order in U they give the same con- 
tribution, =F<? 2 (£//2) 2 /(2i — w p h) 2 , to the integral over the 
spectral function of the self-energy. The sum of the two 
contributions, however, is not zero, but one finds, taking 
into account also terms which involve lower powers of U, 
that it equals — 3g 2 as expected from Eq. I|21|) . There- 
fore, also the sum rule over all frequencies in Eq. (|18|) is 
fulfilled. 

The sum rule for the charge-charge response func- 
tion in Eq. (|25|l applies to finite dopings. For the two- 
site model, this implies the uninteresting filling one, for 
which there is no Coulomb interaction. It is interest- 
ing, however, to study the sum rule over all frequencies 
for the half-filled two-site model. As expected, we find 
that the neglect of vertex corrections, incorrectly, gives 
contributions at lu ~ ±U with weight 1/2, respectively. 
Only when the vertex function from Eq. (|32|l is used in 
Eq. pop. x(<l — TTi^) vanishes to lowest order in t/U. 
This is the correct result from Eq. (|23|l . 
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VII. SUMMARY 



We have derived exact sum rules for the electron- 
phonon contribution S op to the electron self-energy of the 
half-filled Holstcin-Hubbard model in the limit of large 
U. In particular, we consider integrations both over all 
frequencies and over frequencies in the photoemission en- 
ergy range. Comparing results for U = oo and U = 0, 
we find identical sum rules when integrating over all fre- 
quencies but a difference by a factor of four when consid- 
ering frequencies corresponding to photoemission only. 
Using different numerical methods, we find that these 
sum rules are relevant also for systems with intermediate 
values of U w 3D that are typically of interest. These 
sum rules should be useful for testing approximate cal- 
culational schemes. 

We have also used sum rules for studying the impor- 
tance of vertex corrections in a diagrammatic approach to 
properties in the Hubbard-Holstein model. For a weakly 
doped Hubbard-Holstein model in the large-?/ limit, the 
phonon self-energy n is strongly reduced compared to 
the non-interacting case. This is described by a sum 
rule which integrates Im n over a finite frequency range 
such that transitions between the Hubbard bands are not 
involved. This sum rule is satisfied if properly dressed 
Green's functions are used to calculate the phonon self- 
energy, even if vertex corrections are neglected. The en- 
ergy dependence of Im n(q, to) could, nevertheless, be 
substantially influenced by vertex corrections. For the 
half-filled system, we have to integrate Im n over all fre- 
quencies to obtain a nontrivial sum rule. This sum rule 
is only satisfied if vertex corrections are included. 

The sum rule for Im S cp , integrating over frequencies 
corresponding to photoemission only, is violated by a fac- 
tor of four if vertex corrections are neglected. These re- 
sults have been illustrated by an explicit calculation of 
the vertex function in a two-site model. We have studied 
integrated quantities where all values of |q| and u> enter, 
both in terms of their relative ratio and their absolute 
magnitude. Therefore, our findings cannot be directly 
compared with previous ones which focused on the static 
limi1i^4& or on small |q| and ojm^ Our results show that 
the inclusion of vertex corrections can be essential to cor- 
rectly describe effects of electron-phonon interaction in 
strongly correlated systems. 
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APPENDIX A: DERIVATION OF SUM RULES 
USING SPECTRAL MOMENTS 

1. Total sum rule 

To derive the sum rule for the total integrated weight of 
the spectral density of the electron self-energy, Eq. 
we expand [z — lo)^ 1 in Eq. in powers of 1/z. One 
obtains 

oo M ( m ) 

G(M)=£^r, (Al) 
where the moments of the spectral density are defined as 

/oo 
duj u m A(kLu). (A2) 
-oo 

On the other hand, using the Heisenberg equations of 
motion for the time-dependent operators in the defini- 
tion of the spectral density, these moments can also be 
obtained fromi^ 

M^=({C^ct a }), (A3) 

with CO = [O, H}. [O, O'] and {O, O'} denote the com- 
mutator and the anticommutator of two operators O. O' , 
respectively. 

For the Hubbard model with electron-phonon interac- 
tion which was defined in Eqs. CJ-JHJ, one obtains from 
Eq. 123 

Mf - 1, (A4) 
M L 1] = ^ + U{rn- a ) + g^oib, + b\), (A5) 

and 

= 4 + ^ (n ._ ff) + ^| 9q | 2( | 6q + 6 t_ q |2 ) 
q 



+2e k ?/(n l _ CT > + 2e k5q=0 (6 4 + bj) 




where £k, Pqo-, and hi have already been defined after 
Eqs. © and JTSJ. 

When the 1/z expansion of the self-energy, 

oo p(m) 

s ( k ^) = E^r> ( A6 ) 

m=0 

is inserted into Eq. (jHJ , a comparison with Eq. IjAljl leads 
to 

Cf = Mf /Mf - (mW/<>) 2 . (A7) 

As c£ corresponds to the zeroth moment of the spectral 
density of E(k, z), we arrive at the sum rule, Eq. HIGH , 
when Eqs. HA4 |I -I|A6 |I are used in Eq. (|A7jl . 
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2. Sum rules for (I)PES self-energies 

The derivation of sum rules for the (I)PES self-energies 
S ± (k, z) can be done in full analogy with the previ- 
ous section. Only the moments of the (I)PES spectra 
A^(k.,uj) are now obtained from 

M +,(m) = ((£ ™ Ck> t ff)5 (A8) 

M^ m) = (cU m O (A9) 



instead of Eq. I|A3|I . Because Eqs. (|A8|) and (|X9|) do 

not contain an anticommutator, much more complicated 
results are obtained for these moments. We do not list 
them explicitly. It turns out, however, that they simplify 
considerably once the limit U — * oo is taken. Using the 
analog to Eq. (|A7|I and focusing on the electron-phonon 
contribution to the self-energies, one then arrives at the 
result in Eq. fTUf>. 
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